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Abstract 

We study spin dependent phonomena in superconducting junctions in both ballistic and diffusive regimes. For 
ballistic junctions we study both ferromagnet / s- and d-wave superconductor junctions and two dimensional 
electron gas / s-wave superconductor junctions with Rashba spin-orbit coupling. It is shown that the exchange 
field alway suppresses the conductance while the Rashba spin-orbit coupling can enhance it. In the latter part 
of the article we study the diffusive ferromagnet / insulator / s- and d-wave superconductor (DF/I/S) junctions, 
where the proximity effect can be enhanced by the exchange field in contrast to common belief. This resonant 
proximity effect in these junctions is studied for various situations: Conductance of the junction and density of 
states of the DF are calculated by changing the heights of the insulating barriers at the interfaces, the magnitudes 
of the resistance in DF, the exchange field in DF, the transparencies of the insulating barriers and the angle 
between the normal to the interface and the crystal axis of d-wave superconductors a. It is shown that the 
resonant proximity effect originating from the exchange field in DF strongly influences the tunneling conductance 
and density of states. We clarify the followings: for s-wave junctions, a sharp zero bias conductance peak (ZBCP) 
appears due to the resonant proximity effect. The magnitude of this ZBCP can exceed its value in normal states 
in contrast to the one observed in diffusive normal metal / superconductor junctions. We find similar structures 
to the conductance in the density of states. For d-wave junctions at a — 0, we also find a result similar to that in 
s-wave junctions. The magnitude of the resonant ZBCP at a — can exceed the one at a/ir = 0.25 due to the 
formation of the mid gap Andreev resonant states. 
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1. Introduction 

In normal metal / supercunductor (N/S) junctions, a unique scattering process occurs in low energy transport: 
Andreev reflection (AR)[1]. The AR is a process that an electron injected from N with energy below the energy 
gap A is converted into a reflected hole. Taking the AR into account, Blonder, Tinkham and Klapwijk (BTK) 
proposed the formula for the calculation of the tunneling conductance[2]. It revealed the gap like structure or the 
doubling of tunnelilng conductance due to the AR. This method was extended to normal metal / unconventional 
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superconductor (N/US) junctions[3]. It is shown that a zero bias conductance peak (ZBCP) appears when the 
mid gap Andreev resonant state (MARS) is formed due to the anisotropy of US. 

The BTK theory was also extended to ferromagnet / superconductor (F/S) or ferromagnet / unconventional 
superconductor (F/US) junctions[4] and used to estimate the spin polarization of the F layer experimentally 
[5,6,7]. In F/S junctions, AR is suppressed because the retro-reflectivity is broken by the spin-polarization in the 
F layer [8]. To clarify spin dependent transport phenomena is important to fabricate a new device manipulating 
electron's spin. Nowadays, there are many works about charge transport of electrons relevant to electron's spin. 

Among recent works, many efforts have been devoted to study the effect of spin-orbit coupling on transport 
properties of two dimensional electron gas (2DEG)[9,10,11,12]. The pioneering work by Datta and Das suggested 
the way to control the precession of the spins of electrons by the Rashba spin-orbit coupling (RSOC)[13] in 
F/2DEG/F junctions[14]. This spin-orbit coupling depends on the applied field and can be tuned by a gate 
voltage. It also gives the off-diagonal elements of the velocity operator[15]. There are several works about spin 
dependent transport in the presence of RSOC [16,17]. 

The RSOC induces an energy splitting, but the energy splitting doesn't break the time reversal symmetry unlike 
an exchange splitting in ferromagnet. Therefore transport properties in 2DEG/S junctions may be qualitatively 
different from those in F/S junctions. As far as we know, in 2DEG/S junctions the effect of RSOC on transport 
phenomena is not studied well. Recent experimental and theoretical advances in spintronics stimulate us to 
challange this problem. We illustrate the two kind of splittigs in Fig. 1. 

The first purpose of this article is to calculate the tunneling conductance in F/S and 2DEG/S junctions and 
clarify how the exchange field and the RSOC affect it. We think the obtained results are useful for a better 
understanding of related experiments in mesoscopic F/S and 2DEG/S junctions. 

Zeeman vs. Rashba splitting 

Zeeman splitting Rashba splitting 




Figure 1. Zeeman vs. Rashba splitting. 

On the other hand, in diffusive junctions the physics is clearly different from that in ballistic junctions. In 
diffusive normal metal / superconductor (DN/S) junctions, proximity effect plays an important role in the low 
energy transport. The phase coherence between incoming electrons and Andreev reflected holes persists in DN 
at a mesoscopic length scale and results in strong interference effects on the probability of AR [18]. One of the 
striking experimental manifestations is the zero bias conductance peak (ZBCP) [19,20,21,22,23,24,25,26,27,28,29]. 

A quasiclassical Green's function theory [30, 31, 32] is often applied to the charge transport in DN/S junctions. 
Volkov, Zaitsev and Klapwijk (VZK) solved the Usadel equations [33], and showed that this ZBCP is due to the 
enhancement of the pair amplitude in DN by the proximity effect [34]. VZK applied the Kupriyanov and Lukichev 
(KL) boundary condition for the Keldysh-Nambu Green's function [35]. Stimulated by the VZK theory, several 
authors studied the charge transport in various junctions [36,37,38,39,40,41,42,43,44]. 

Recently one of the authors [45] developed the VZK theory for s-wave superconductors using more general 
boundary conditions provided by the circuit theory of Nazarov [46]. The boundary conditions coincide with 
the KL conditions when a connector is diffusive with low transparent coefficients, while the BTK theory [2] is 
reproduced in the ballistic regime. The extended VZK theory [45] produced a crossover from a ZBCP to a zero bias 
conductance dip (ZBCD). These phenomena are relevant for the actual junctions in which the barrier transparency 
is not necessarily small. 
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The formation of the MARS at the interface of unconventional superconductors [47,3,48,49] also generates the 
ZBCP as mentioned above. The generalized VZK theory was recently applied to unconventional superconducting 
junctions [50,51,52]. The formation of the MARS is naturally taken into account in this approach. It was demon- 
strated that the formation of MARS in DN/d-wave superconductor (DN/D) junctions strongly competes with the 
proximity effect. 

Above theories treat spin independent phenomena in diffusive junctions. Calculations of tunneling conductance 
in the presence of the magnetic impurities in DN /S junctions were performed in Ref[53,54,34]. Spin dependent 
transport is also realized in ferromagnet / superconductor junctions. 

In diffusive ferromagnet / superconductor (DF/S) junctions Cooper pairs penetrate into the DF layer from the 
S layer and have a nonzero momentum by the exchange field[55,56,57]. This property produces many interesting 
phenomena[58,59,60,61,62,63,64,65,66,67,68,69,70,71,72]. One interesting consequence of the oscillations of the 
pair amplitude is the spatially damped oscillating behavior of the density of states (DOS) in a ferromagnet 
predicted theoretically [73,74,75,76]. In the ferromagnet the exchange field usually breaks the induced Cooper 
pairs. But for a weak exchange field the pair amplitude can be enhanced and the energy dependent DOS can have 
a zero-energy peak[77]. The DOS has been studied extensively [74, 77, 78, 79] but the condition for the appearance 
of the DOS peak was not studied systematically. We studied the conditions for the appearance of such anomaly, 
i.e., strong enhancement of the proximity effect and found two conditions corresponding to weak proximity effect 
and strong one[80]. Since DOS is a fundamental quantity, this resonant proximity effect can influence various 
transport phenomena. 

Another purpose of the present article is to study the influence of the resonant proximity effect by the exchange 
field on the tunneling conductance and the DOS in DF/ s- and d-wave superconductor junctions with Nazarov's 
boundary conditions. Weak exchange field is realized in recent experiments with, e.g., Ni doped Pd[79] or magnetic 
semiconductor. Thus our results may be observed in experiments. In the latter part of this article we calculate the 
tunneling conductance and the density of states in normal metal / insulator / diffusive ferromagnet / insulator / 
s- and d-wave superconductor (N/I/DF/I/S) junctions for various parameters such as the heights of the insulating 
barriers at the interfaces, resistance Rd in DF, the exchange field h in DF, the Thouless energy Eth in DF, the 
transparencies of the insulating barriers and the angle between the normal to the interface and the crystal axis of 
d-wave superconductors a. Throughout the paper we confine ourselves to zero temperature. 

The organization of this paper is as follows. In sections II and III, we will provide the detailed derivation of 
the expression for the normalized tunneling conductance and the results of calculations are presented for various 
types of junctions, in ballistic and diffusive junctions respectively. In section IV, the summary of the obtained 
results is given. 



2. Ballistic junctions 

We consider F/S and F/d-wave superconductor (F/D) junctions. We use the same method as in Ref. [4] and the 
same notations. In the following ] (J.) denotes majority (minority) spin. The F/US interface located at x = (the 
j/-axis) has an infinitely narrow insulating barrier described by the delta function U(x) = H8(x). As a model of 
the ferromagnet we apply the Stoner model with the exchange potential U. The pair potential matrix we consider 
is given by 



where 6 denotes the direction of motions of quasiparticles measured from the normal to the interface. Below 
we consider s- and d-wave superconductors. The pair potentials are given by Ajf = — Afj = A for s-wave 
superconductors and A|j = — Afj. = Acos[2(# — a)] for d-wave superconductors where a denotes the angle 
between the normal to the interface and the crystal axis of d -wave superconductors. Here A denotes the energy 
gap. 

Applying the BTK method[2,4], we obtain the conductance ffsT(l) f° r U P (down) spin quasiparticle represented 
in the form: 




(1) 



0~St = a N1 



l- ir+r-i^i-q^ + q^ir+i 2 



?9{0c-\0\) 



|1 - r + r_^/l - ajvtV 1 - °"ivi exp [i (ip l - ipf)]\ 



3 



i - ir + r_ 



1 - T+T-y/l - a N f exp [i (ipx - </j t ) 



;Q(\e\-6 c ) 



i — |r + r_ | (l - o-jv T ) + q-jv T |r + 

l - r+r-V 1 - o-jvtV 1 = °ni exp [i (v3| - <pi) 



rQ(e c -\0\) 



(2) 
(3) 



exp(ipj) = 



1 — A_ + iZ e 



s/1 - ctni (1 + A- 



oxp(-iv9 T ) 



1 - A+ - iZ e 



r± 



^/l-a^ (l + A+ + iZ e ) 
A± (0) 



i?+^i5 2 -|A ± (0)| 2 ' 



(4) 
(5) 



Z ° - ^7e> Z - Wk 



1,1 * and 9c = cos 1 J with quasiparticle energy E, effective mass m, Fermi wavenumber 



kF and Fermi energy Ef- In the above 6>(:r) is the Heaviside step function and <Jni(i) denotes the conductance 
for up (down) spin quasiparticle in the normal state: 
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with 



Normalized conductance is expressed as 



ctt — 



S-jl ddcos, 9 (asf + asi) 
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(7) 



Next we consider a ballistic 2DEG/S junctions. The 2DEG/S interface located at x = (along the y-axis) has 
an infinitely narrow insulating barrier described by the delta function U(x) = HS(x). The effective Hamiltonian 
with RSOC is given by 



H = 



( Cfc iAfc_e(-a:) A9(j) \ 
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(8) 



with fc± = k x ± Cfc — ^7 ~ ^f)i Fermi wave number fop, Rashba coupling constant A. 
Velocity operator in the x-direction is given by[15] 
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Wave function ip(x) for x < (2DEG region) is represented using eigenfunctions of the Hamiltonian: 
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for an injection wave with wave number fci(2) where fci 



+ 



k 2 



+ 



f#) +k% and 



fci(2)± = fci(2)e * 1(2) - ii(2) and bi( 2 ) are AR coefficients. Ci( 2 ) and di( 2 ) are normal reflection coefficients. 9i(2) is 
an angle of the wave with wave number fci( 2 ) with respect to the interface normal. 

Similarly for x > ip(x) (S region) is given by the linear combination of the eigenfunctions. Note that since 
the translational symmetry holds for the y-direction, the momenta parallel to the interface are conserved: k y = 
kF sin# = fci sin#i = k 2 sin #2- 

The wave function follows the boundary conditions[15]: 



v x ip (x)\ x=+0 - v x tp {x)\ x= _ = 



h 2mU 
mi K 2 
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Applying BTK theory to our calculation, we obtain the dimensionless conductance represented in the form: 



cr s = Al 



J \ {K-a. + \ ai \ 2 K 21 + |6i| 2 K 12 \ 2 i - |ci| 2 K 21 - |di| 2 K 12 \ 2 i] cosOdt 
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(13) 



with Ki2 = l + fe, if: 
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(14) 



Ni and A 2 are normalized density of states for wave number fci and fc 2 respectively. The critical angle 0c is 
defined as cos <9c = y fr^. 

(j n is given by the conductance for normal states, i.e., as for A = 0. We define normalized conductance as 
err = crs/ un and a parameter /3 as f3 = 2 ^^ . 

First we study the difference between the effect of the Zeeman splitting and that of Rashba splitting. We plot 
the tunneling condutance for superconducting states, as for F/S junctions in (a)-(c) and for 2DEG/S junctions 
in (d)-(f) of Fig. 2 with Z = 10 in (a) and (d), Z = 1 in (b) and (e), and Z = in (c) and (f). The exchange 
field suppresses as independently of Z as shown in (a)-(c). This is because the AR probability is reduced by the 
exchange field. On the other hand the dependence of as on j3 at zero voltage is qualitatively different. In (d)-(f) 
we show the dependence of as on j3 at zero voltage for various Z. For Z = 10 it has an exponential dependence on 
(5 but its magnitude is very small. For Z — 1 it has a reentrant behavior as a function of (5. For Z = it decreases 
linearly as a function of /3. The AR probability is enhanced by the RSOC at Z = 10. 
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Next we will study the F/D junctions. The normalized tunneling conduntace gt as a function of bias voltage 
V is plotted in Fig. 3 for a/ir = 0.25 and various exchange field with Z = 10 in (a) and Z = in (b). For 
Z — 10 a ZBCP appears due to the formation of the MARS as shown in (a). As the exchange field increases, <tt 
is suppressed. Similar plots at Z = are shown in (b). We can find that <tt decreases with the increase of the 
exchange field. 
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Figure 2. Tunneling condutance for superconducting states at zero voltage as a function of the exchange field in F/S junctions (left 
panels) and RSOC in 2DEG/S junctions (right panels) with Z — 10 in (a) and (d), Z — 1 in (b) and (c), and Z — in (c) and (f). 
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Figure 3. Normalized tunneling conductance in F/D junctions at a/-K — 0.25 with Z — 10 in (a) and Z — in (b). 



3. Diffusive junctions 



We consider a junction consisting of normal and superconducting reservoirs connected by a quasi-one-dimensional 
diffusive ferromagnet conductor (DF) with a length L much larger than the mean free path. The interface between 
the DF conductor and the S electrode has a resistance Rb while the DF/N interface has a resistance R' b . The 
positions of the DF/N interface and the DF/S interface are denoted as x = and x — L, respectively. We model 
infinitely narrow insulating barriers by the delta function U(x) = H8(x — L) +H'8(x). The resulting transparency 
of the junctions T m and T' m are given by T m = 4cos 2 0/(4cos 2 <j) + Z 2 ) and T' m = 4cos 2 0/(4cos 2 <f> + Z' 2 ), where 
Z = 2H/vf and Z' = 2H' /vf are dimensionless constants and (j> is the injection angle measured from the interface 
normal to the junction and vf is Fermi velocity. 
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We apply the quasiclassical Keldysh formalism in the following calculation of the tunneling conductance. The 
4x4 Green's functions in N, DF and S are denoted by Gq(x), Gi(x) and ($2 (a?) respectively where the Keldysh 
component Ko,i,z{x) is given by Ki(x) — Ri(x)fi(x) — fi(x)Ai(x) with retarded component Ri(x), advanced 
component Ai(x) — —R*(x) using distribution function fi(x)(i = 0,1,2). In the above, Ro(x) is expressed by 
Ro(x) = T3 and /o(a?) = fio + Tsfto- R-2(x) is expressed by R,2{x) = g% + /fa with g — e/\/e 2 — A 2 and 
/ = A/VA 2 - e 2 , where T2 and T3 are the Pauli matrices, and s denotes the quasiparticle energy measured from 
the Fermi energy and fo(x) = tanh[e/(2T)] in thermal equilibrium with temperature T. We put the electrical 
potential zero in the S-electrode. In this case the spatial dependence of Gi (x) in DF is determined by the static 
Usadel equation [33], 

D ^[Gi(x)^ 1 ]+i[H,G 1 (x)] = (15) 
with the diffusion constant D in DF. Here H is given by 
f Ho \ 

H=\ . , (16) 
\ Ho J 

with H = (e — (+)ft)f 3 for majority(minority) spin where h denotes the exchange field. Note that we assume 
a weak ferromagnet and neglect the difference of Fermi velocity between majority spin and minority spin. The 
Nazarov's generalized boundary condition for Gi(x) at the DF/S interface is given by Refs.[45,51]. We also use 
Nazarov's generalized boundary condition for Gi(x) at the DF/N interface: 

L ( ^ dG! _ .-1 , 2T4[Go(0-),G 1 (0+)] 

Wd ( Gl — )^ 0+ --R b < B >' B - 4 + T4([ g o( o_), Gi(0+ )] + -2)- 

The average over the various angles of injected particles at the interface is defined as 

d<b cos d>B(<b) 

< b{4>) > (,) = J -;' 2 v — — — (is) 

/:^ 2 #TW(^)cos0 
with B(cf>) = B and T^(<f>) = T$ . The resistance of the interface R b is given by 



J#> = 4° — = . (19) 

/:^ 2 ^TW(0)cos0 



Here 7?g'' is Sharvin resistance, which in three-dimensional case is given by Rq ' 1 = e 2 k 2 F S^ /^-K 2 ), where k F 
is the Fermi wave-vector and Sc ^ is the constriction area. Note that the area Sc'' is in general not equal to the 
crossection area Sd of the normal conductor, therefore Si'' / Sa is independent parameter of our theory. This allows 
to vary Rd/Rb independently of T m . In real physical situation, the assumption Si'^ < Sd means that only a part of 
the actual flat DN/S interface (having area Si'^) is conducting whether it a single conducting region or a series of 
such regions. These conducting regions are not constrictions in a standard sense - we don't assume the narrowing 
of the total cross-section, but rather that only the part of the cross-section is conducting. 
The electric current per one spin is expressed using G\{x) as 

L f d eTriM Gl (x)^±n (20) 



8eRd J dx 


where (Gi(x) dG g^ ) K denotes the Keldysh component of (Gi(x) dG g^ )■ In the actual calculation it is convenient 
to use the standard ^-parameterization where the retarded Green's fucntion Ri(x) is expressed as Ri(x) = 
fa cos 9(x) + T2 sin 9(x). The parameter 8(x) is a measure of the proximity effect in DF. 

The distribution function /1 (x) is given by /1 (x) = fi(x) + %ft{x). In the above, ft(x) is the relevant distribution 
function which determines the conductance of the junction we are now concentrating on. From the retarded or 
advanced component of the Usadel equation, the spatial dependence of 6(x) is determined by the following equation 



d 2 

D q-2 6{ - x ) + 2i ( e - sm[6>(z)] = (21) 
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for majority(minority) spin, while for the Keldysh component we obtain 

D^-[^lco S h 2 e im (x)} = 0. (22) 
At x = 0, since fto is the distribution function in the normal electrode, it is satisfied with 

fto = ^{tanh[(e + eV)/{2T)] - tanh[(e - eV)/(2T)]}. (23) 

Next we focus on the boundary condition at the DF/N interface. Taking the retarded part of Eq. (17), we obtain 

L d0(x) , < F >' 2T m sm9 

■-^ = 7^ — ^TTT^ „„„« ' ( M > 



R d dx i% (2 - T^)+T4 cos < 

with 6> = 0(0+). 

On the other hand, from the Keldysh part of Eq. (17), we obtain 

L ( df t , 2 < Jm >' (fto - /t(0+)) . T^A; + 27^(2 - r^)Real{cos go} 
_(_ )cosh * im(a!) | x=0+= _ >7kl = _______ (25) 

Ai = (1+ | cos6> | 2 + | sin(9 | 2 ). (26) 
After some calculations we obtain the following final result for the current 



el 
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oo 

^ f de f - _ (27) 



<I b0 > T L JO cosh^ e irn (x) T </(,!>' 

Then the differential resistance i? per one spin projection at zero temperature is given by 



2% , •_>//., / : '/.) 2./?; 



with 



_£ / da: + _* (28) 
< J w > ' L ./ cosh 2 6» lm (x) <J 6 i>' 
o 

, T 2 A! + 2r m (2-r m )A 2 

60 2 I (2 - T m ) + T m [gcos6>L + /sin6»z,] | 2 ' 1 yj 

Ai = (1+ | cos8 L | 2 + | sinfe | 2 )(| g | 2 + | / | 2 +1) + 4Imag[/ff*]Imag[cos0i, sin6»2], (30) 

A 2 = Real{g(cos6> L + cos6>£) +f(sin# L +sin6»£)}. (31) 

This is an extended version of the VZK formula[34]. In the above 9i m (x) and 9l denote the imaginary part of 
9(x) and 9(L-) respectively. Then the total tunneling conductance in the superconducting state as(eV) is given 
by as(eV) = _t J. V-^- The local normalized DOS N(e) in the DF layer is given by 

N(e) = i^Rccos6»(a;). (32) 
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It is important to note that in the present approach, according to the circuit theory, R d /R b can be varied 
independently of Tm , i.e., independently of Z^'\ since one can change the magnitude of the constriction area 
Sc^ independently. In other words, Rd/R^ is no more proportional to TaJ(L/l), where T^y is the averaged 
transmissivity of the barrier and I is the mean free path in the diffusive region. Based on this fact, we can choose 
Rd/R^ and as independent parameters. 

In the following, we will discuss the normalized tunneling conductance or(eV) = as(eV) / 'crjv(eV) where (Tjv(eV) 
is the tunneling conductance in the normal state given by <Tjv(eV) = a at = l/(-R<2 + Rb + R't)- 

Now we study the influence of the resonant proximity effect on tunneling conductance as well as the DOS in the 
DF region. The resonant proximity effect was discussed in Ref. [80] and can be characterized as follows. When the 
proximity effect is weak (Rd/Rb <S 1), the resonant condition is given by Rd/Rb ~ 2h/ETh due to the exchange 
splitting of DOS in different spin subbands. When the proximity effect is strong (R d /Rb 3> 1), the condition is 
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given by Eth ~ h and is realized when the length of a ferromagnet is equal to the coherence length £f = V ' D Jh. 
We choose Rd/Rb — 1 as a typical value to study the weak proximity regime. We also choose Rd/Rb = 5 to study 
the strong one. We fix Z — Z' = 3 because these parameter don't change the results qualitatively and consider 
the case of high barrier at the N/DF interface, Rd/R' b ~ 0.1, in order to enhance the proximity effect. 

Let us first choose the weak proximity regime and relatively small Thouless energy, Eth I A = 0.01. In this case 
the resonant condition is satisfied for h/Exh = 0.5. In Fig. 4 we show the tunneling conductance for Rd/Rb = 1, 
ETh/A = 0.01 and various H/Eth in (a). The ZBCD occur due to the proximity effect for h = 0. For h/ Eth = 0.5, 
the resonant ZBCP appears and split into two peaks or dips at eV ~ ±ft with increasing h/Exh- The value of 
the resonant ZBCP exceeds unity. Note that ZBCP due to the conventional proximity effect in DN/S junctions 
is always less than unity [21,45] and therefore is qualitatively different from the resonant ZBCP in the DF/S 
junctions. 

The corresponding normalized DOS of the DF is shown in (b) and (c) of Fig. 4. Note that in the DN/S junctions, 
the proximity effect is almost independent on Z parameter [45]. We have checked numerically that this also holds 
for the proximity effect in DF/S junctions. Figure 4 displays the DOS for Z — 3, Rd/Rb = 1 and Et%/^ = 0.01 
with (b) h/Exh = and (c) h/Exh = 0.5 corresponding to the resonant condition. For h — 0, a sharp dip appears 
at zero energy over the whole DF region. For nonzero energy, the DOS is almost unity and spatially independent. 
For h/Exh = 0.5 a zero energy peak appears in the region of DF near the DF/N interface. This peak is responsible 
for the large ZBCP shown in (a). Therefore, the ZBCP in DF/S junctions has different physical origin compared 
to the one in DN/S junctions. 




Figure 4. Normalized tunneling conductance (a) and the DOS ((b) and (c)) with Rd/Rb — 1 and Eth/& — 0.01. (b) K/Eth — 
and (c) h/E Th = 0.5. 



Next we choose the strong proximity regime and relatively small Thouless energy, _Bt/i/A = 0.01. In the present 
case, the resonant ZBCP is expected for h/Erh = 1. Figure 5 displays the tunneling conductance for Rd/Rb = 5 
and J^Th/A = 0.01 and various h/Exh in (a). We can find resonant ZBCP and splitting of the peak as shown in 
(a). The corresponding DOS is shown in (b) K/Eth = and (c) K/Eth = 1. For h = 0, a sharp dip appears at 
zero energy. For finite energy the DOS is almost unity and spatially independent. For K/Eth — 1 a peak occurs 
at zero energy in the range of x near the DF/N interface. We can find a similar structure in the corresponding 
conductance as shown in Fig. 5(a). The DOS around zero energy is strongly suppressed at the DF/S interface 
[x — L) compared to the one in Fig. 4. 

Let us proceed to study the d-wave junctions both for weak and strong proximity regimes. In this case, depend- 
ing on the orientation angle a, the proximity effect is drastically changed: As a increases the proximity effect 
reduces[50,51]. For a = we can expect similar results to those in the s- wave junctions since proximity effect exists 
while the MARS is absent. In contrast, the tunneling conductance for large a is almost independent of H/Eth- 
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Figure 5. Normalized tunneling conductance (a) and the DOS ((b) an d ( c )) with R^/R^ — 5 and Eth/^- — 0.01. (b) h/Ep^ — 
and (c) H/Eth — 1- 

Especially, the conductance is independent of h for a/ir = 0.25 due to the complete absence of the proximity 
effect. There exist two different origins for ZBCP in DF/D junctions: the ZBCP by resonant proximity effect 
peculiar to DF and the ZBCP by the MARS formed at DF/D interface. When a deviates from 0, the MARS are 
formed at the interface. At the same time, the proximity effect is suppressed due to the competition between the 
proximity effect and the MARS. Therefore the MARS provide the dominant contribution to the ZBCP compared 
to the resonant proximity effect in DF, as will be discussed below. 

First we choose the weak proximity regime where the resonant condition is Ji/Eth = 0.5. Figure 6 displays the 
tunneling conductance for Rd/Rb = 1, -Et/i/A = 0.01 and various a with (a) h/Em = and (b) h/Exh = 0.5. 
For h = 0, ZBCD appears for a/n = due to the proximity effect as in the case of the s-wave junctions while 
ZBCP appears for a/n = 0.25 due to the formation of the MARS (Fig. 6 (a)). For h/E Th = 0.5, the height of 
the ZBCP by the resonant proximity effect for a — exceeds the one by MARS for a/n = 0.25 (Fig. 6 (b)) in 
contrast to the ballistic junctions where the ZBCP for a/n = 0.25 is most enhanced[3]. 

We also study the DOS of the DF for the same parameters as those of Fig. 6(b) as shown in (c) a/n — 
and (d) a/n = 0.125 in Fig. 6. The line shapes of the LDOS near the DF/S interface are qualitatively similar 
to the tunneling conductance. For a/n — a zero-energy peak appears as in the case of s-wave junctions. With 
increasing a, the DOS around the zero energy is suppressed due to the reduction of the proximity effect. The 
extreme case is a/n — 0.25, where the DOS is always unity since the proximity effect is completely absent. 

Next we look at the junctions for the strong proximity regime. In Fig. 7 we show the tunneling conductance for 
Rd/Rb = 5, Eth/A = 0.01 and various a with (a) h/ETh = and (b) h/Exh = 1- In this case we also find the 
ZBCP for a — by the resonant proximity effect. The height of the ZBCP is suppressed as a increases as shown 
in Figs. 7(b). 

The corresponding DOS of the DF to the case of (b) in Fig. 7 is shown in Fig. 7 with (c) a/n — and (d) 
a/n = 0.125. For a — the structure of the DOS at x — reflects the line shapes of the tunneling conductance. 
With increasing a the zero energy peak of the DOS becomes suppressed. The DOS at the DF/S interface (x = L) 
are drastically suppressed compared to those in Fig. 6. 
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Figure 7. Normalized tunneling conductance ((a) and (b)) and the DOS ((c) and (d)) for d-wavc superconductors with R&jR 
and E Th /A = 0.01. (a) h/E Th = 0, (b) h/E Th = 1 (c) h/E Th = 1 and q/tt = and (d) h/E Th = 1 and q/tt = 0.125. 



11 



4. Conclusions 



In this article we studied the tunneling conductance in F/S, F/D and 2DEG/S junctions in ballistic regime. We 
extended the BTK formula and calculated the tunneling conductance of the junctions. We clarified the following 
point: 

1. The exchange field always suppresses the conductance in F/S and F/D juntions. 

2. In 2DEG/S jucntions, for low insulating barrier the RSOC suppresses the tunneling conductance while for 
high insulating barrier it can slightly enhance the tunneling conductance. We also found a reentrant behavior of 
the conductance at zero voltage as a function of RSOC for intermediate insulating barrier strength. The results 
give the possibility to control the AR probability by a gate voltage. We believe that the obtained results are useful 
for a better understanding of related experiments of mesoscopic F/S and 2DEG/S junctions. 

In the latter of the present paper, a detailed theoretical study of the tunneling conductance and the density of 
states in normal metal / diffusive ferromagnet / s- and d-wave superconductor junctions is presented. We clarified 
that the resonant proximity effect strongly influences the tunneling conductance and the density of states. There 
are several points which have been clarified in this paper: 

3. For s-wave junctions, due to the resonant proximity effect, a sharp ZBCP appears for small Eth ■ We showed 
that the mechanism of the ZBCP in DF/S junctions is essentially different from that in DN/S junctions and is 
due to the strong enhancement of DOS at a certain value of the exchange field. As a result, the magnitude of the 
ZBCP in DF/S junctions can exceed unity in contrast to that in DN/S junctions. 

4. For d-wave junctions at a — 0, similar to the s-wave case, the sharp ZBCP is formed when the resonant 
condition is satisfied. At finite misorientation angle a, the MARS contribute to the conductance when Rd/Rb <C 1 
and Z > 1. With the increase of a the contribution of the resonant proximity effect becomes smaller while the 
MARS dominate the conductance. As a result, for sufficiently large a ZBCP exists independently of whether the 
resonant condition is satisfied or not. In the opposite case of the weak barrier, Rd/Rb 3> 1, the contribution of 
MARS is negligible and ZBCP appears only when the resonant condition is satisfied. 

An interesting problem is a calculation of the tunneling conductance in normal metal / diffusive ferromagnet / 
p-wave superconductor junctions because interesting phenomena were predicted in diffusive normal metal / p-wave 
superconductor junctions[52]. We will perform it in the near future. 
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